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ABSTRACT 

Dark energy has a dramatic effect on the dynamics of the universe, causing the recently 
discovered acceleration of the expansion. The dynamics are also central to the behav- 
ior of the growth of large scale structure, offering the possibility that observations 
of structure formation provide a sensitive probe of the cosmology and dark energy 
characteristics. In particular, dark energy with a time varying equation of state can 
have an influence on structure formation stretching back well into the matter domi- 
nated epoch. We analyze this impact, first calculating the linear perturbation results, 
including those for weak gravitational lensing. These dynamical models possess defi- 
nite observable differences from constant equation of state models. Then we present a 
large scale numerical simulation of structure formation, including the largest volume 
to date involving a time varying equation of state. We find the halo mass function 
is well described by the Jenkins et al mass function formula. We also show how to 
interpret modifications of the Friedmann equation in terms of a time variable equation 
of state. The results presented here provide steps toward realistic computation of the 
effect of dark energy in cosmological probes involving large scale structure, such as 
cluster counts, Sunyaev-Zel'dovich effect, or weak gravitational lensing. 
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1 INTRODUCTION 



Direct dynamical measurement of the expansion of the 
universe through the Type la supernova distance-redshift 
method discovered th a t the expansion is accelerating 
JPerlmutter et alJ fl999l: iRiess et all I1998T) . This has wide 
reaching implications with respect to the fate of the uni- 
verse, its dominant constituent, and the nature of funda- 
mental physics. Some 70% of the total energy density acts 
like a dark energy with strongly negative pressure. Sub- 
sequent observations of the cosmic microwave background 
(CMB) power spect rum and of large s cale structure give a 
concordant picture feond et alJl2002t iPercival et all 120021: 
ISpereel et atll2003ft . 

Mapping the expansion history of the universe offers a 
way to gain insights into the mysterious dark energy. For 
example, characterizing its equation of state (pressure to 
energy density ratio) behavior gives insight on the prop- 
erties of the high energy physics scalar field potential. Dis- 
tance measures, notably the supernova method, have proved 
adept at beginning to constrain the energy density and equa- 
tion of state of the dark energy. Recent limits, combined 
with CMB or large scale structure information, within the 
low redshift approximation of a constant equation of state 
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(EOS) ratio w, gi ve —1.61 < w < —0.78 at 95% confidence 
llKnop et al.ll2003) . Great improvements should occur in the 
next dec ade with, e.g., the Su pernova/ Acceleration Probe 
(SNAP; lAldering et al] 2002) dedicated dark energy pro- 
gram, in particular extending to constraints on the generi- 
cally expected time varying function w(z). 

Several other cosmological probes look promising, 
though their systematic uncertainties are less well defined. 
But the eventual synergy of independent and complemen- 
tary methods should prove powerful in revealing the nature 
of dark energy. Some of these probes offer the opportunity 
to measure fairly directly the expansion rate behavior H(z), 
rather than just the distance which involves this quantity 
through a redshift integral. In turn, H(z) involves an in- 
tegral over the equation of state w(z). One cannot however 
naively assume that a limited measure of H(z) is better than 
knowing the distance over a wide redshift range. This was 
demon strated f or the c osmic shear, or Alcock-Paczyhski, ef- 
fect in fLindei <l2002b[) and the baryon osci llation probe in 
lLinderll2003bh fsee lHuterer fc Turnerll200ll for a variety of 
other methods). But such information can prove valuable in 
complement with precision distance data. 

Furthermore, methods involving the growth of large 
scale structure appear, initially at least, to possess sensi- 
tivity to the cosmic equation of state. This enters through 
the actual growth of density fluctuations, i.e. the balance 
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of attractive gravitational stability with the dynamic fric- 
tion of the expansion, through its evolution with redshift or 
time as the dynamics changes under the influence of dark 
energy, and through the cosmic volume available in which 
structures form. However little rigorous work has been done 
on the full impact of dark energy other than for the time 
independent cosmological constant model. This especially 
applies to the virtually universal time varying EOS models, 
where not just the energy density but the function w(z) is 
time dependent. Indeed a constant w model (other than the 
cosmological constant w = —1) functions only as a crude 
approximation, unsuitable for next generation data that ex- 
tends beyond z m 0.5 or that seeks to combine complemen- 
tary probes. Moreover, the important and revelatory physics 
responsible for the accelerating universe appears in the field 
dynamics - the time variation w' ~ dw/dz. 

If we desire to take advantage of the power of large 
scale structure formation as a probe of dark energy, we 
must include sufficient realism in the model that we form 
a consistent picture of the underlying physics. This includes 
both systematic uncertainties in the astrophysics and ob- 
servations, and a treatment of time variation in the dark 
energy EOS. Structure based methods such as weak lensing, 
Sunyaev-Zel'dovich distance measures, and cluster counts all 
require knowledge of how structure formation behaves in the 
presence of realistic dark energy. In Section|5|we present the 
key role of w' in describing the dynamics of the expansion 
and in Section [3] solve the growth equation in linear pertur- 
bation theory. Section 2] describes the numerical simulation 
of large scale structure in the presence of dark energy and 
Section discusses the results. Implications and an outline 
of future research are presented in Section 



EOS and the time variation w' = dw/dln(l+z)\ z =i = w a /2, 
allows consideration of data extending to z > 1 and presents 
an excellent approximation to slow roll scalar field dark en- 
ergy models (Linder 2002a). Thus we can use the observa- 
tions in a well defined manner to investigate the fundamental 
physics manifesting in the EOS time variation. 

However, a time varying equation of state is more gen- 
eral than this. Some theories have been proposed that ex- 
plain the acceleration not through a scalar field dark energy 
but through modifications of the Friedmann equations them- 
selves by alternative theories of gravitation, e.g. arising from 
extra dimensions, or by highly speculative components such 
as the Chap lygin ga s or qua ntum or higher dimensional cor- 
rections. In iLindeil i2003al) . a formalism was presented to 
use supernova distance data to constrain these scenarios by 
means of mapping the expansion history a(t) directly, rather 
than using an intermediate parametrisation w(z). While this 
is valid, and probably physically preferred, we note here an 
alternate interpretation. 

Consider Eq. ©. The dark energy term really just de- 
scribes our ignorance concerning the physical mechanism 
leading to the observed effects of acceleration in the expan- 
sion, i.e. an increase in the expansion rate. Let us instead 
write this as 



H 2 {z)/Hl = Q, m {l + z) 



6H 2 /H 2 , 



(3) 



where now we encapsulate any modification to the Fried- 
mann equation of general relativity in the last term. That 
is, we take a very empirical approach: all we have observed 
for sure is a certain energy density due to matter, Q m , and 
consequences of the expansion rate H(z). 

We can now write the deceleration parameter generally 



2 w' IS EVERYWHERE 

The expansion rate of the universe, H(z) — a/a, where a(t) 
is the scale factor, enters into both the kinematics and dy- 
namics of the cosmological model. Distances are integrals of 
the proper or conformal time; for example in a flat universe 
(assumed throughout) the comoving distance is 



r(z) = / dt/a(t) = / dz/H(z) 



(1) 



where the redshift z = a -1 — 1. The angular diameter dis- 
tance is just r a = r/(l + z) and the luminosity distance 
is ri = (1 + z) r. The expansion rate, or Hubble parame- 
ter, H(z) also enters into dynamical quantities such as the 
growth of structure through a "Hubble drag" term. 

Within the flat universe, dark energy picture, the Fried- 
mann equations give the expansion rate as 

rr2, \/tj2 /-, . \3 . /■, o \ 3 /„ dln(l + z')[l + w(z')] 

H (z)/H = U m (l+z) +(1— ilm)e J o . 



(2) 



where Ho is the Hubble constant, the present value of the 
Hubble parameter, Q m is the dimensionless matter density 
today (so the dark energy density is 1 — Q m ), and w(z) 
is the generically time dependent dark energy equation of 
state. While each model of dark energy has a particular form 
for w(z), in order to compare models on e usuall y adopt s 
a parametrisation. The one proposed by iLindeJ (l2003af) : 
w(z) — wq + w a z/(l + z), with wo the present value of the 
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If we interpret the modified expansion rate as being due to a 
w(z) as appearing in Eq. - whether or not it has anything 
to do with a scalar field - then we find 

1 dlnSH 2 



^^ 1+ 3rfffi(TTI)- (5) 

This now defines an effective, time v arying equation of state 
(something similar was noticed bv lAlam etafl 12003ft . Of 
course it reduces to the usual result in the scalar field dark 
energy case. But this goes to illustrate the centrality of a 
time varying w(z) - or something that looks just like it - 
for probing cosmological models. To fix the EOS to be a 
constant w rather than a varying function w(z) is highly 
non-generic, an unjustified assumption and a frequently poor 
approximation, and can blind us to important physics. 



3 LINEAR PERTURBATION THEORY 

The growth of structure depends sensitively on the expan- 
sion rate of the universe. For example, the solution to the 
classic Jeans instability in a static space shows exponen- 
tial growth under gravity, while this gets reduced to only a 
power law behavior in time in an expanding space-time. The 
perturbations are sourced in the gravitational instability of 
slightly denser regions having correspondingly greater grav- 
itational attractions and thus further increasing in density; 
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this is opposed by an effective friction, or "Hubble drag" 
term, due to the expansion. 

Since dark energy affects the expansion rate, one ex- 
pects to see an influence due to the friction. That this can be 
substantial comes from our experience that open universes, 
which effectively have a component with EOS w — —1/3, 
can shut off the growth of structure when the curvature en- 
ergy density dominates over the matter density. Similar re- 
sults hold for a cosmological constant dominated universe. 
Thus for general dark energy models the state of structure 
formation at various redshifts could probe the equation of 
state, and even its time variation (we fix the sound speed 
to its canonical value, c s = 1). We first present a general, 
pedagogical analysis of the influence of EOS on growth in 
the linear perturbation regime, then specifically calculate 
the results for various dark energy models. 



3.1 Growth of Linear Perturbations 

On scales smaller than the horizon, t he dark energy 
component is expected to be sm ooth jMa et alJ Il999t 
iDave. Caldwell, fc SteinhardtteOOST) so we only consider per- 
turbations to the matter. Then the growth equation becomes 



5" 



5 + 2H5 ■ 
- q)a~ 1 5' - 



(3/2)Q m a' 2 5 



(6) 
(7) 



where 8 is the fractional matter density perturbation, q is the 
deceleration parameter, dot denotes a time derivative and 
prime a derivative with respect to scale factor a. One can 
readily see that growth in a universe with Qtot and fi_i/3 
in a component characterized by w = —1/3 behaves like 
growth in a universe with no such component but Q' to t — 
Qtot — Q-i/3- So a flat universe with aw = —1/3 component 
acts like an open universe. 

We can write the growth equation in terms of the gen- 
eral EOS w(z), which as we saw in Sj5] can also represent 
modifications to the framework of the theory. Because the 
interpretation of the source and drag terms is so straight- 
forward, this generalization broadly (but not always: see the 
object lesson within standard gravitation in the Appendix) 
carries over to the growth dynamics. Defining the growth as 
the ratio of the perturbation amplitude at some scale fac- 
tor relative to som e initial scale factor, D = S(a) /5(a,i), the 
equation becomes (lLinderlll988T) 



1 - 



w(a) 



1 + X(a) 



_ 3 X(a) 
a 2 1 + X(a) a 2 ~ ' 



-3 j" din a' w(a' ) 



(8) 

(9) 
(10) 



where Eq. Ogives the general case of X and Ea. 1101 puts it 
in terms of the time dependent scalar field equation of state. 

The variable X is the ratio of the matter density to the 
dark energy density, and the growth equation holds even if 
flatness does not, so we do not set X/(l + X) — Q m (a). As 
X gets small the source term vanishes and growth cannot 
generally proceed. Note that this is the reason why (sub- 
horizon dark matter) perturbations cannot effectively grow 
in the radiation dominated epoch, not from any overriding 



influence of the Hubble drag term. (Of course we have not in- 
cluded coupling between components; radiation-baryon cou- 
pling prevents gravitational instability in the baryon per- 
turbations). From Eq. JSJ, one sees that the friction term 
opposing the growth is proportional to 1 — w when matter 
is not dominant. So in the radiation epoch the drag is less 
than in an open or accelerating epoch. It is only for those 
cases when w < 1/3, as comes from solving the characteris- 
tic equation, that the growth is shut off by the friction term. 
(Indeed growth can occur for w > 1/3.) 

One can readily verify that for large X one recovers the 
matter dominated behavior D ~ a. It is convenient to re- 
move this trend, common to all the considered cosmological 
models at high redshift, and define a variable G = D/a. The 
evolution equation for this "normalized" growth is 



G" + 



3 w(a) 



2 1 + X(a) 



G' 3 1- w(a) G _ 
a 2 l+X(a) a? 1 ' 



Through the cosmological Poisson equation G is related to 
both the gravitational potential and peculiar velocity fields. 

This formalism readily allows incorporation of dark en- 
ergy with a time varying EOS, or the alternate equation of 
motion Eq. Figure illustrates the solutions for a vari- 
ety of dark energy models. Note that models with constant 
equation of state have modest differences from the cosmo- 
logical constant case, ~ 7% for w — —0.8. This is still larger 
than the differences caused by changing the matter density 
by ±0.02, shown by the dotted curves flanking the cosmo- 
logical constant model. 

Dark energy whose EOS is negatively evolving, w' < 0, 
e.g. acts more like a cosmological constant in the past but 
with a less negative EOS today, is almost indistinguishable 
in its linear growth predictions from the cosmological con- 
stant. This is because the dark energy contribution to the 
total energy density only becomes significant at late times 
in these models. This property holds as well for dark energy 
models with EOS w < — 1. 

However, the class of models with w' > exhibits dra- 
matically different behavior. This includes what we chose as 
our fiducial model to test tim e varying EOS, the sup ergrav- 
ity inspired model SUGRA of lBrax fc Martini lll999l) that is 
well fit by wq = - 0.82, w a = 0.58 (i.e. w' = +0.29). Note 
that Linder (2002a) has shown that the fit is good to 3% in 
w at z = 1.7 and to 0.2% in distance to the last scattering 
surface at z — 1089. 

Models with w' > show definite differences in the 
growth of structure from a cosmological constant universe. 
For a model like SUGRA (here actually wo = —0.8, w' = 
0.3) the growth disparity reaches 26% at the present day. 
The figure shows two other points of interest. Comparing 
SUGRA to the constant EOS model that matches its dis- 
tance to the last scattering surface, and so effectively mimics 
the time varying model as far as CMB data is concerned (ex- 
cept at low multipoles), the difference in total growth is 12%. 
That is, the time varying model and the constant, effective 
model can be distinguished through large scale structure in- 
formation, despite being largely degenerate in their CMB 
power spectra. This is a promising sign. 

The other interesting element is that even at high red- 
shift the time varying model with w' > possesses dif- 
ferent behavior from the cosmological constant and other 
constant w models. This shows the influence of dark en- 
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Figure 1. Growth factor of linear density perturbations plotted 
vs. scale factor a = (1 + z)~ 1 for six dark energy models with 
fim = 0.3. Dotted curves show the effect of changing O m by 
±0.02 for the cosmological constant case. Curves labeled w e ff 
use the constant equation of state that gives the same distance 
to the last scattering surface as the bracketed models, and hence 
are degenerate with respect to the CMB. 



ergy at early times is very different from the cosmological 
constant, which quickly becomes dynamically negligible for 
2 > 1. Such a characteristic offers the possibility that evi- 
dence for time varying EOS dark energy might be found as 
well in the low multipole region of the CMB power spectrum 
fsee lCaldwell et alJ l'2003. for a discussion of the influence of 
"early quintessence" on the CMB through the integrated 
Sachs- Wolfe effect). 

So three avenues appear to be open for the detection 
of physically important properties of dark energy from its 
influence on large scale structure: 1) linear growth rate, 2) 
nonlinear structure formation and evolution, and 3) large 
angle power in the CMB power spectrum. In the remainder 
of this section we investigate further the first avenue and 
then proceed in !jl]to discuss numerical simulations of the 
second possibility fsee lBenabed fc Bernadeau Il200ll for an 
an alytic atte mpt). The third approach has been addressed 
by_]Hu| (2002) with not very optimi stic conclusions (but see 
ICoorav. Huterer fc Baumannll2003T) . 

Somewhat different results for the linear growth rate 
appear if we change the measurement technique. Suppose 
that rather than normalizing the density perturbations by 
their high redshift behavior (currently corresponding to 
COBE/WMAP normalisation of the matter power spec- 
trum, though Planck will provide great improvement), we 
calibrate them by their present amplitude. This is like fix- 
ing as, the power on the interface of linear/nonlinear scales. 
While observations have not yet determined this precisely, 
we can explore the consequences. As seen in Fig. now 
the models are difficult to distinguish at low redshifts and 
there is near degeneracy in the growth factor between the 



Figure 2. Same as Fig.^but with the growth factors normalized 
by their values today. The present growth, related to the mass 
variance erg, is not yet precisely known, however. 

time varying model and the effective constant w model that 
matches it with respect to the CMB. 

The issue of the optimal way to measure the growth fac- 
tor through the matter power spectrum is an area of ongoing 
research, increasingly important with the future advanced 
large scale structure surveys. Perhaps the most realistic ob- 
servable for now is the ratio of growth factors at different 
redshifts, a measure of the evolution of structure. This will 
carry correspondingly less content since it is only a relative 
measure, without information from the absolute level. Such 
a ratio could be read off from either Fig. Qor|2]and indeed 
can be seen to vary little with cosmological model. For ex- 
ample, the evolution between z = 2 and z — 1 in the models 
with (wq, Wo) = (—0.8, 0) or (—0.8, 0.3) agrees with the evo- 
lution in the cosmological constant model to 2% or 5%. This 
casts strong doubt on the idea that the growth evolution of 
large scale structure by itself (without a precise and robust 
high redshift, CMB normalisation) is useful as a probe of 
dark energy. 

3.2 Sensitivity of Linear Growth Rate to Dark 
Energy 

To understand the use of the evolution of linear matter den- 
sity perturbations for probing the dark energy, we need to 
study not just the gross differences in the curves in Fig. 
but the details of how they depend on dark energy proper- 
ties and other cosmologica l parameters. We use the Fisher 
matrix formalism (see, e.g.. lTegmark et~ Il998h to plot the 
sensitivity of the growth to the parameters f2 ra , wo , and w a 
in Fig. |!|] The sensitivity increases toward low redshift since 
this corresponds to more time for the differing growth dy- 
namics to take effect. 

Besides the sensitivity, the degeneracy between param- 
eters is a crucial aspect to the usefulness of the probe. The 
similarity of the shapes of the wo and w a curves guarantees 
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Figure 3. Logarithmic sensitivity of the growth factor to the 
cosmological parameters, perturbed around the cosmological con- 
stant model, as a function of redshift. Such a plot is useful in find- 
ing lower bounds to parameter estimation errors and in indicating 
degeneracies between parameters. 



a strong degeneracy between them, apart from further inter- 
action with the value of f2 m . Combined with the relatively 
low accuracy achievable on the growth factor from observa- 
tions, this leads to the growth factor by itself - even with 
the initial amplitude known - not serving as a precision cos- 
mological probe. For example, even if all other parameters 
were fixed, a 5% determination of S/a at z — 1, say, would 
only constrain wo to ±0.35. And degeneracies strongly am- 
plify the uncertainty. Additionally, the degeneracy direction 
in the wo — w a plane is roughly aligned with the contours 
from the CMB, so the growth factor possesses little comple- 
mentarity with the CMB and only a modest amount with 
the supernova distance measure. 

For the simulated data sets we consider a measurement 
of the growth factor to 5% (equivalent to 10% determina- 
tion of the linear power spectrum) at three values of the scale 
factor (equivalent to z = 0.67, 1.5, 4), SNAP supernova mea- 
surements (including systematic uncertainties), and Planck 
CMB determination of the angular distance to the last scat- 
tering surface (normalisation of the primordial power is im- 
plicit in the growth factor measurement). The results, shown 
in Fig. |3] are not sensitive to the exact redshifts chosen for 
the growth factor estimates. Addition of growth information 
improves constraint of wo, w a by less than 2%. This holds 
as well if the growth factor is normalized to z = 0. 

However, as was found for the baryon oscillation method 
l)Linderl2003lJ) . the situation changes dramatically when the 
fiducial dark energy model is not taken to be the cosmologi- 
cal constant but the SUGRA model with time varying EOS. 
Now the growth factor under the influence of the dynamic 
field has a different cosmological dependence and possesses 
substantial complementarity with both the supernova and 
CMB data, as shown in Fig. [£] The parameter uncertainties 
using all three methods in synergy become a(wa) = 0.04, 



Figure 4. Parameter estimations (68% confidence level) of the 
present equation of state and its time variation, marginalizing 
over the matter density Q m with a prior of <r(Q m ) = 0.03. Dot- 
ted curves use the growth factor normalized to its value today. 
The growth factor alone (not shown) is poor in constraining the 
cosmological model, and offers little complementary leverage for 
the proposed SNAP supernova survey, and none in addition to 
SNAP plus the Planck CMB survey. 



a(w') = 0.05 (recall that w' = w a /2). That gives a 6a de- 
tection of time variation in the EOS! 

A further interesting characteristic of the growth factor 
is a sufficiently distinct cosmological dynamics dependence 
to break some degeneracies outside the dark energy frame- 
work. For example, the brane world model w ith crossover 
scale Hor c = 1.43 discussed in iLinderl <2003alh difficult to 
distinguish from a w = —0.7 dark energy model through the 
distance relation (< 0.5% difference), does differ by ~ 4% 
in the growth factor (though not when normalized at low 
redshift). This illustrates the utility of the SH 2 formalism of 
321 for testing the cosmological framework. 

3.3 Growth Rate and Weak Gravitational Lensing 

One of the main applications of the growth rate is to cosmo- 
logical observables that statistically characterize large scale 
structure. In [21 we address how the ingredient of the linear 
growth factor enters description of nonlinear structure. But 
here we make a brief, illustrative foray into the linear regime 
of the weak lensing shear used to map out the matter dis- 
tribution (including dark matter) through its gravitational 
deflection of light from distant sources. 

Observations of such shear are becoming increasingly 
useful cosmological tools and depend on the primordial mat- 
ter density power spectrum, the growth rate, and geomet- 
ric distance factors. We examine the simplified, though cen- 
tral, quantity of the shear linear le nsing power spectrum (cf. 
I Jain fc Selial3ll99'itlHuteredl2002h : 



G, 



dzi 
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Figure 5. Same as Fig. ||] but for dark energy following the 
SUGRA model, with present equation of state too = —0.82 and 
time variation w a = 0.58. In this case, information on the linear 
growth factor adds valuable constraints. 



Figure 6. Parameter estimations (68% confidence level) of dark 
energy properties using future 5% measurement of the weak grav- 
itational lensing linear power spectrum in addition to SNAP su- 
pernova distances and Planck CMB data. The lensing power has 
good complementarity with supernovae, though little with the 
CMB. Thus supernovae plus lensing should provide a valuable 
crosscheck on supernovae plus CMB results. 



Here Zi is the redshift of the lensing mass, z s of the source 
galaxy, r s the comoving distance to the source, and r; s the 
comoving distance between the source and the lens. For sim- 
plicity we used the linear matter power spectrum ~ k G 2 
and fix z 3 . The latter corresponds to having tomographic 
information. 

We find that this lensing combination alone cannot 
place useful constraints on the dark energy but it does have 
good complementarity with the supernova distance measure- 
ments. Such lensing data strongly improves the determina- 
tion of fl m , as traditionally expected, and reduces uncer- 
tainty in the time variation w a slightly more dramatically 
than CMB information. However it has little complemen- 
tarity with CMB data, and so the limits are not much fur- 
ther improved (1%) on inclusion of both lensing and CMB 
data. Fig.|S|shows the case of 5% determination of the weak 
lensing shear power spectrum at z s — 1.5. Results are not 
very sensitive to the exact redshift, or combination, chosen. 
The SUGRA case leads to similar conclusions, though there 
the weak lensing has more complementarity with the CMB 
and adding the lensing information tightens the constraints 
on u>o and w a from supernovae and CMB by 34%. Thus 
the SNAP weak lensing program in combination with the 
SNAP supernova survey provides an important crosscheck 
on the results from SNAP supernovae plus CMB, plus the 
possibility of additional improvements. 

Of course the probe may be further strengthened by us- 
ing information from the nonlinear part of the power spec- 
trum. To employ the nonlinear regime, one needs to carry 
out numerical simulations including a time varying EOS or 
use fitting formulas valid for such cases, neither of which 
previously existed. In the remainder of the paper we present 
and discuss such a simulation. 



4 A SIMULATION OF LARGE SCALE 
STRUCTURE FORMATION 

So far in this paper we have concentrated on the dependence 
of the linear growth factor on the nature of the dark energy. 
In this section we turn our attention to the influence of the 
dark energy on the growth of non-linear large-scale structure 
- specifically the dark matter halo mass function. This mea- 
sure of the number density of objects as a function of mass 
is important in a wide range of large-scale structure areas, 
including the use as a probe of the cosmological model. 

Theoretically the dark matter halo mass function has 
been explor ed in CDM models main l y through n-body sim - 
ulations e.g.lEfstathiou et alJ (Il98l) :jLacev fc Cole I <1994): 
iGross et alJ iGovernato et al J <1999l) : I Jenkins et alJ 

<200ll) : IWhitd (l200ll. l2002t) . Nearly all simulations to date 
focussing on this issue have modelled a universe in which the 
dark energy is static - a cosmolological constant A - and the 
matter component is exclusively dark - and have ignored 
the fluid and dissipative nature of the baryonic component. 
The baryons are taken into account when computing the 
input matter transfer function used in making the initial 
conditions for the simulation. In the simulation itself it is 
assumed that the baryons behave as a pressureless compo- 
nent which mirrors the dark matter precisely. Modelling of 
the formation of individual haloes where a baryonic com- 
ponent is included, for example the S anta Barbara cluster 
comparison project jFrenk et al Jl999t) . lends support to the 
idea that one can determine the mass function for virialised 
dark matter haloes well without including the baryons at 
least for the high mass end. 

At the same time as progress has been made studying 
the mass function using n-body simulations, analytic formu- 
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lae have been developed to describe the mass functions. The 
most influential work was pioneered by IPress fc Schechteil 
which utilises the spherical top-hat model of col- 
lapse. Comparisons of the Press-Schechter mass function 
to n-body simulations by a number of authors have noted 
that the Press-Schechter mass function tends to overesti- 
mate the mass function at low masses and underestimate 
the high mass end. This discrepancy has led to the develop- 
ment of im proved analytic mass fun c tion formulae, such as 
the w o rk of ISheth fc Tormenl Jl999h ; ISheth. Mo fc Tormenl 
j200ll) : ISheth fc Tormenl J2002I) fS-T hereafter), based on a 
model which accounts for ellipsoidal rather than spherical 
collapse. However in this paper we will conc entrate on the 
' unive rsal' mass function formula presented in I Jenkins et alJ 
i200ll) (J01 hereafter), which is an empirical fit to a wide 
range of simulation data (of which those used by S-T form 
a subset). As we will show in the next section this for- 
mula does actually predict the mass functions in a SUGRA 
quintessence model without needing modification. 

Both simulations and the analytic formulae show for 
CDM models with Gaussian power spectra that the high 
mass end of the halo mass function is very sensitive to the 
precise normalisation of the power spectrum. If one can by 
some means determine the halo mass function over a range of 
redshifts then in principle one can find how the linear growth 
factor evolves as a function of redshift for the universe. 

The dark matter halo mass function is not directly ob- 
servable, though future observations exploiting weak grav- 
itational lensing will make progress in this direction. How- 
ever it is believed that the largest dark matter haloes are 
recognisable because they host galaxy clusters. Galaxy clus- 
ters contain, as well as galaxies, copious amounts of x-ray 
emitting gas. The presence of x-ray emitting gas not only 
signals the presence of the cluster halo which helps to find 
and count them, but can also be exploited to estimate the 
actual mass of the surrounding dark matter halo in a va- 
riety of ways. Estimates of the halo mass function using 
observational data of galaxy clusters have been made and 
utilised by many authors to esti mate the value of the matter 
power quantity as loca l ly - e .g. [White, Efstathiou & Frenk 
jl993Myiana fc Liddlel l)l996h . lEke. Cole fc Frenk I 1)199*6^ 
Seliak (2002). Observations of clusters over a range of red- 
shifts potentially enables the evolution of the growth fac- 
tor to be determined. Two difficulties to surmount in using 
clusters as cosmological probes are the relation of the ob- 
servables to theoretical characteristics, e.g. mass, and then 
the extraction of cosmological parameters free from astro- 
physical factors and degeneracies. The typical uncertainty 
in the theoretically determined mass functions are around 
10%. However making accurate mass estimates of galaxy 
clusters is far from easy and comparing an observationally 
determined mass function with the relatively clean theoreti- 
cal mass function requires extensive and detailed modelling 
- it is not our aim in this paper to discuss these important 
aspects. 

As mentioned above the mass function formula given 
in J01 is a fit to the mass functions from a range of CDM 
simulations. These simulations include models with a cos- 
mological constant (w = — 1) and Open models which have 
an effective value of w — —1/3. The values of interest for w 
for the dark energy are typically within the range that has 
already been studied so it would not seem too surprising if 



the mass function for models with constant values of w inter- 
mediate between —1/3 and —1 are accurately predicted by 
existing mass function formulae. However it is interesting to 
ask how the mass function of a model where the value of w is 
rapidly changi ng compares to th e predictions. The SUGRA 
model of iBrax fc Martini l)l999T) is an attractive model to 
choose for this reason because the value of w evolves rea- 
sonably rapidly. 

The part of the mass function which is most sensitive to 
the linear growth factor is the high mass end where the mass 
function falls off steeply. One would similarly expect that it 
would be this region of the mass function that would show 
the greatest sensitivity to a time variation in w. We therefore 
have designed our simulation to model a large volume of 
space so that we can measure the abundance of very massive, 
but rare, clusters. 

As far as we are aware the only simulations with 
a quint essence component with a changing EOS are 
those of iKlypin et~aT1 J2003t) who also m odel the SUGRA 
quinte ssence together with a model by iRatra fc Peeblesl 
(1988), which has a much more gently changing value of 
w, and some models with constant w. Their simulations 
have very good spatial and mass resolution, allowing them 
to study the properties of individual dark matter haloes in 
some detail. However their simulation cubes are relatively 
small compared to ours: even their largest simulation cube 
is some 66 times smaller in volume than our cube. The re- 
sults the authors report for the mass functions determined 
from their simulations appear to be fully consistent with 
our own results. Our simulations are designed to give better 
statistics at the high mass end of the mass function and are 
in this sense complementary. 

Below we describe the n-body code, the parameter 
choices for the simulation and the initial conditions for the 
simulation. The results themselves are reported in the sec- 
tion |S] 

4.1 Code details. 

We have used the publ icly available parallel code GAD- 
GET (vl.l) described in lSpringel. Yoshida fc White! (1200 ll) 
to perform the n-body simulation. We modified the code 
to allow the inclusion of a quintessence component with an 
EOS of the form w(a) — wo + (1 — a)w a , where a is the 
expansion factor defined so that a = 1 at the present epoch. 
The only modifications required to the code were to up- 
date expressions for the Hubble parameter to include the 
quintessence component. We verified these modifications by 
checking that the modified code can reproduce the correct 
linear growth rates on a series of test simulations. We will 
only report the results of this test for the simulation pre- 
sented here. 

We determined the linear growth rates from the n- 
body simulation by measuring the power associated with 
the longest wavelength modes of the simulation box for the 
initial conditions and for each output. The longest wave- 
length modes are least affected by non-linear effects. For 
comparison we also computed the growth rates by numerical 
integration of equation lilt . Between z=38.9 and z=2 the 
agreement in the growth factors measured from the simula- 
tion and computed by numerical quadrature is better than 
0.2%. This is a non-trivial test that the code can handle a 
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quintessence component because the SUGRA model has an 
EOS w — —0.25 at the upper redshift and by redshift 2 the 
quintessence component already has a significant dynami- 
cal influence: for example the value of G(a) is 0.84 rather 
than unity as it would be in a matter dominated universe. 
The long wavelength modes grow slightly more slowly than 
the linear theory prediction, due we b elieve to non-linear 
effects fsee lBaugh fe Efstathiou1 [T999'). By z=Q, when the 
long wavelength density fluctuations themselves have an am- 
plitude of a few percent, the agreement is better than 0.8%. 
By comparison this difference is 25 times smaller than the 
difference of the growth rates of the quintessence model and 
ACDM over the same interval. 

4.2 Simulation details. 

We have designed the simulation with two aims in mind: 
(i) to determine the high end of the halo mass function, 
and (ii) to match as closely as possible, at redshift zero, 
the parameters of an existing ACDM simulation so that 
a direct comparison can be made. The most suitable ref- 
erence ACDM simulation for determining the high mass 
end of the halo mass function is th e ACD M Hubble vol- 
ume run described in lEvrard et al. I (120021) which models 
a cube, 3000/i _1 Mpc on a side. We have therefore selected 
the same values for the matter density, Q m , the particle 
mass and the gravitational softening length as the ACDM 
Hubble simulation. We also match very closely the power 
spectrum shape and amplitude of the linearly extrapolated 
power spectrum at redshift zero to those of the ACDM 
Hubble simulation. Our simulation volume is significantly 
smaller than the ACDM Hubble simulation at 648fe _1 Mpc 
on a side. 

We have chosen a quintessence component with a vari- 
able value of w. As discussed earlier, the choice wo = —0.82, 
w a = 0.58, gives a good fit to the SUGRA quintessence 
model. We will refer to our simulation as the SUGRA- 
QCDM model. The parameters of the simulation are listed 
in table^ 

The initial conditions were generated by the serial ver- 
sion of the code that was used to generate t he initial con- 
dition s for the Hubble volume simulations dEvrard et al. I 
2002). The initial conditions are created from an initially 
uniform particle distribution (in this case a glass distribu- 
tion generated in the way described bv lWhitelll996l) by per- 
turbing the particles to give the desir ed power spectru m and 
using the Zel'dovich approximation jZerdovichlll970|) to as- 
sign velocities to each particle which are proportional to the 
displacements. The code needed to be altered to give the 
correct constant of proportionality, which depends on the 
logarithmic growth dlnD/dlna, when assigning the veloc- 
ities. The difference in this constant between an Einstein- 
deSitter (Q m = 1) and quintessence model at the start red- 
shift is small - only ~ 1.5% - though this is much larger than 
the difference between a matter dominated and cosmological 
constant model. 



5 SIMULATION RESULTS 

The simulation was started at z=38.9 and outputs were 
made at redshifts z = 3, 2, 1.5, 1, 0.5, 0.25 and 0. 



Table 1. Parameters of the SUGRA-QCDM n-body simulation - 
the mass resolution, softening and power spectrum normalisation 
are chosen to match the ACDM Hubble volume simulation de- 
scribed in Evrard et al 2002. <rg is the RMS of the linear density 
field smoothed with a top-hat filter of 8fo —1 Mpc in radius. 



Parameter Value 



n m 0.30 

w -0.82 

w a 0.58 

<t 8 (z = 0) 0.90 

Box size//i _1 Mpc 648 

Number of particles 10077696 

Particle mass/10 12 h- 1 M Q 2.25 

Comoving softening lcngth/h _1 Mpc 0.1 



5.1 The mass function at z=0 

To begin with we will compare the mass functions at redshift 
zero. The z=0 linear power spectrum of the SUGRA-QCDM 
simulation has the same shape and amplitude as the ACDM 
Hubble simulation by d esign. Analytical fit t ing fo rmulae for 
mass functions such as IPress fc Schechtej (Il974h . S-T and 
Ijenkins et al.l i200lft . all predict that the CDM mass func- 
tion depends primarily on the linear power spectrum. We 
would therefore expect if one extrapolates these results to 
quintessence models that our two models should have very 
similar mass functions at redshift zero. 

One potential complication in comparing two different 
cosmological models lies in the problem of how to define the 
haloes in a consistent way so that they can be compared. 
A common approach to this problem is to use the spherical 
top-hat collapse model to give guidance as t o the expected 
collapse overdensity of virial is ed ob jects - e.g. lLacev fc Cole I 
(1993ft : lEke. Cole fc Frenk I dl996ft . Following these argu- 
ments would lead to different values for the overdensity for 
haloes selected in ACDM and SUGRA-QCDM. 

However we can circumvent these issues by adopting 
the halo definition used in J01. T he haloes are defin ed using 
the friends-of-friends algorithm JPavis et alJll985ft with a 
linking length of b=0.2. J01 showed that using this way of 
defining haloes it is possible to fit the mass functions for 
CDM models with a wide range of cosmological parameters 
and redshifts with a single 'universal' fitting formula which 
is accurate to better than 20%. It is this simplicity which 
motivates our choice of halo definition. In this subsection we 
will compare the SUGRA-QCDM and ACDM simulations at 
z=0 directly. 

Our ability to determine the mass function is limited by 
a number of considerations. At the high mass end, where the 
objects themselves are resolved with the largest numbers of 
particles, the main limitation in determining the mass func- 
tion is the volume of the box. We plot the mass functions 
only up to the point where the Poisson error reaches 10%. 
At the low mass end the situation is less simple. One would 
expect determination of the mass function to be less and less 
reliable as one resolves haloes with fewer and fewer particles. 
We take a limit of 20 particles as the minimum number of 
particles. This limit is supported by tests presented in ap- 
pendix A of JO 1 where they compare the mass functions from 
simulations with the same mass resolution as our SUGRA- 
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Friends-of-friends mass function (b = 0.2) 
z = 0. 



ACDM Hubble 



\M,z) = 
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log 10 (M/h-'M s ) 

Figure 7. Top panel: the dark matter halo mass function defined 
using the fricnds-of-fricnds algorithm with a linking length of 0.2. 
The mass functions of ACDM Hubble simulations and SUGRA- 
QCDM that at z=0 have nearly identical power spectra are very 
similar. The mass functions are plotted for haloes of 20 particles 
and more and up to an upper mass cut-off where the Poisson 
uncertainty in the mass function first exceeds 10%. Bottom panel: 
the solid curve shows the difference between the SUGRA-QCDM 
mass function and the ACDM Hubble volume mass function. The 
differences are below 10%. 



QCDM simulation to those from simulations with superior 
mass resolution. 

In figure □ we show the SUGRA-QCDM and ACDM 
mass function evaluated with the FOF(b=0.2) group finder. 
The mass functions from both simulations have been 
smoothed using an identical Gaussian with an RMS width of 
0.08 dex. The mass functions agree very well - the differences 
are below 10%. Thi s is an encouraging result ( a similar re- 
sult is also seen in iKlvpin et al"l [2003) and suggests that it 
should be possible to predict the mass functions for this and 
other quintessence models with a smoothly varying value 
of w with a precision of ~ 10% from a knowledge of just 
the linearly evolved power spectrum. In the next section we 
will compare the SUGRA-QCDM mass function with the 
universal mass function formula of J01 directly for several 
redshifts. 



5.2 The mass function at all redshifts 

Here we will plot the combined mass functions at redshifts 
3, 2, 1.5, 1, 0.5, 0.25 and for the SUGRA-QCDM run and 
compare them with the J01 mass function formula. But we 
first repeat a few definitions, adapted from J01, which are 
required further on. 

It is convenient to use as an effective mass variable, 
Inc , where a(M) is the RMS of the linear density field 
smoothed with a top-hat filter containing mass M at the 
mean density. This is defined as: 



k 2 P(k)W 2 (k- M)dk, 



(13) 



where D(z) is the growth factor of linear perturbations nor- 
malised so that D(0) = 1, P(k) is the linear power spectrum 
at redshift zero and W(k;M) is the Fourier-space represen- 
tation of a real-space top-hat filter enclosing mass M at the 
mean density of the universe. 

We define the mass function /(ln<r -1 ,z) through: 



/(In a 



-l _ M dn(M,z) 
, z) = 



po diner" 



(14) 



where n(M, z) is the abundance of halos with mass less than 
M at redshift z, and Po{z) is the mean density of the uni- 
verse at that time. With these definitions we can plot the 
mass functions for any CDM model and at any redshift con- 
veniently onto the In / — lna -1 plane. 

Figure |H1 shows the mass functions for the SUGRA- 
QCDM simulation, plotted in the In / — In a^ 1 plane, for 
redshifts z= 0, 0.25, 0.5, TO, 1.5, 2 and 3. The mass func- 
tions are plotted as dashed lines for haloes of 20 particles 
and more (as in J01) and cut-off at the high mass end when 
the Poisson errors reach 10%. The simulation curves have 
been smoothed using a Gaussian with an RMS of 0.05 dex. 
The solid line shows the J01 universal mass function fit given 
by: 



/(In a" 1 ) =0.315 exp [ 



lna -1 +0.61I 3 - 8 ! 



(15) 



valid over the range —1.2 ^ lna -1 ^ 1.05, and the flanking 
dotted lines denote a 20% uncertainty about this fit. These 
curves are similarly convolved with a 0.05 dex Gaussian. 

Note that the universal nature does not imply that 
different models predict the same observable consequences, 
only that they can be treated by the same parametrization. 
Each model is distinguished by its particular a 2 (M, z) rela- 
tion. 

The main result of this section is that the SUGRA- 
QCDM mass function is well fit by the J01 formula - at least 
to the same degree of accuracy as was found for a host of 
other CDM cosmological models. The fit is good despite the 
fact that this model has a changing value of w. It therefore 
seems likely that the J01 should work well for a wide range of 
dark energy models. Although we have used the J01 formula 
in figure El w e could instead have used the formula of S-T. 
IKlvpin et alJ i2003T) find that the S-T formula provides an 
excellent fit to their mass functions. We do find however that 
for lno --1 above about 0.75 that J01 formula is definitely a 
better fit, with the S-T formula tending to overpredict the 
number of haloes. 

It is apparent from figure |H| that the most significant 
deviation away from J01 formula occurs for lner" 1 > 0.85, 
where only the z = 3 output contributes. A closer analysis 
of this curve and other curves shows that the plotted mass 
functions tend to drop away relative to the J01 curve for 
small numbers of particles per group. This drop off is par- 
ticularly apparent for the z = 3 and 2 outputs where the 
mass function is steep. The cause of this feature is unclear 
but is likely numerical in origin. The influence of the dark 
energy component diminishes as one goes to higher redshift 
and the J01 formula works well for matter dominated mod- 
els - so if anything one would expect the degree of fit to 
improve with increasing redshift contrary to what is seen. 
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Figure 8. The halo mass functions for the SUGRA-QCDM simu- 
lation plotted in the / — lnc - 1 plane for redshifts 0,0.25,0.5,1,2,3. 
The smooth solid curve flanked by dotted curves marks the J01 
universal mass function formula. The dotted lines mark a ±20% 
uncertainty in the mass predicted mass function. The SUGRA- 
QCDM mass functions are plotted for groups of 20 or more 
particles and up to the point where the Poisson error becomes 
10%. The SUGRA-QCDM mass function is well fit by the Jenk- 
ins et al formula. Some departure is apparent at high values of 
lncr - 1 > 0.85, but as discussed in the text the origin of this de- 
parture is likely due to numerical limitations rather than being 
real. 

This divergence is much less pronounced at lower redshift 
and is not apparent in the resolution tests presented in ap- 
pendix A of J01 for the Hubble volume simulations which 
are for redshift zero. 

To conclude, firstly we have found in this section that 
the z — mass functions of the Hubble ACDM simulation 
and our own SUGRA-QCDM simulation match to a preci- 
sion of better than 10% in abundance for haloes in the mass 
range of 4.5 x 10 13 -2 x 10 15 h _1 M Q . By design, the z = lin- 
ear power spectrum of the two models was closely matched. 
The results indicate that the mass function depends primar- 
ily on the linear power spectrum and is only very weakly if 
at all dependent on the details of the expansion history. Sec- 
ondly, from an analysis of the SUGRA-QCDM mass func- 
tion at a range of redshifts we find that mass function for- 
mulae such as J01 and S-T provide good estimates of the 
halo mass function for flat quintessence models. The J01 
formula gives a better fit at the high mass end than S-T for 
lna -1 > 0.75. The goodness of fit is comparable with that 
found when comparing these formulae to the mass functions 
determined from the following cosmological models: matter 
dominated, open, and flat models with a cosmological con- 
stant. We are able to verify the accuracy of the fit over range 
-0.2 < In a" 1 < 0.85 from our simulation. 

These results imply that the main observational dis- 
criminators for large-scale structure between cosmological 
models are: the present linear power spectrum, to be fixed 
by wide field surveys; the linear growth factor, discussed 



in the first part of this paper and probed by future deep 
surveys; and distances and volumes measured by expansion 
history mappers such as SNAP. 



6 CONCLUSIONS 

Cosmic structure formation and evolution provides an ad- 
ditional path to exploring the cosmological model besides 
mapping the expansion history. Many of the observational 
probes employing structure rely on the growth dynamics of 
density fluctuations, in either the linear or nonlinear regime. 

Linear perturbations offer the best hope of structure 
measurements with relatively uncomplicated physics and 
clean observations free of many astrophysical entanglements. 
Analyzing the linear growth behavior, we find little lever- 
age on the nature of the dark energy - unless it possesses a 
time varying equation of state. In this latter case the growth 
shows good complementarity with supernova and CMB mea- 
surements. 

We have emphasized that such time variation appears 
generically in proposed dark energy models other than the 
cosmological constant and can also be used to treat modifi- 
cations of the cosmological framework. 

The growth also enters as a contribution to the gravi- 
tational lensing power spectrum, and weak lensing as a cos- 
mological probe promises interesting complementarity and 
crosschecks. Such a method can take advantage of both the 
linear and nonlinear density scales to attempt to balance 
sensitivity with systematic uncertainties. 

For the nonlinear realm of structure, investigations at 
the level of precision necessary require numerical simula- 
tions. We present one of the first, and currently the largest 
volume, N-body simulation with a time varying equation of 
state. As a first application of the data, we calculate the 
halo counts as a function of mass, a quantity relevant to 
forthcoming structure surveys. We find that this is well fit 
by the previous Jenkins mass formula, extending its uni- 
versality to time varying equations of state, at least at the 
20% level. At redshift zero the dynamical dark energy model 
shows agreement to better than 10% with the cosmological 
constant simulation that matches the linear power spectrum 
today, suggesting this serves generally as a central quantity 
in describing structure formation. 

Much work remains for the future to bring observations 
of structure to the level of theoretical and systematic uncer- 
tainty necessary for precision probing of dark energy. But 
the results here lay a foundation to build upon for using not 
only the expansion history of the universe but the growth 
history of structure in a rigorous quest to understand the 
physics behind the accelerating universe. 
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APPENDIX. APPROXIMATION OF THE 
LINEAR GROWTH FACTOR 

Within the formalism of the Bir khoff's theorem argument 
presented bv lPeebiesl Jl98lll993h for the evolution of linear 
density perturbations, one can write a closed form expres- 
sion for the growth. However this does not hold for general 
cosmological models. While it was sufficient and indeed pre- 
scient for the main models considered in 1980, it neglects a 
term arising for a general pressure component (as Peebles 
alludes to); this term is proportional to (1 + 3w))(l + w)Q w 
and so we see that serendipitously the closed form is exact 
for a pure matter model (SCDM, £l w — 0), an open model 
(OCDM, id = —1/3), and a cosmological constant model 
(ACDM, w — —1). But a second order differential equation 
does not generally possess such a quadrature and instead an 
equation like Eq. must be solved. 

Fig. [5] shows the difference between the exact solution 
for the growth factor and the approximation given by the 
closed form 

5 H f a 3 
Gapprox = — fi m — / da(aH) . (16) 

z a Jo 

The differences for the time varying SUGRA model can be 
~ 15%; this propagates into the power spectrum as the 
square, and then further as a possibly substantial bias on 
the cosmological parameters. The situation is better if future 
detailed observations allow use of the growth factor normal- 
ized to the present (~ 1 % difference for SUGRA, t hough 
~ 5% for w = -1.2). See lWang fe Steinhardtl lll99Sft for a 
fitting function involving an integral over the equation of 
state, valid at high redshifts, for slow variations. 
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Figure 9. Comparison of the exact growth factor (solid curves) 
with the closed form approximation, Eq. lltjl fdashed curves). The 
solid curves are labeled according to the vertical ordering, and 
the subscripts indicate the ordering of the dashed curves for those 
models. The approximation is insufficient for precision cosmology. 



